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Abstract 

The harmonic oscillator with dissipation is studied within the framework of the Lindblad 
theory for open quantum systems. By using the Wang-Uhlenbeck method, the Fokker- Planck 
equation, obtained from the master equation for the density operator, is solved for the Wigner 
distribution function, subject to either the Gaussian type or the ^-function type of initial 
conditions. The obtained Wigner functions are two-dimensional Gaussians with different 
widths. Then a closed expression for the density operator is extracted. The entropy of 
the system is subsequently calculated and its temporal behaviour shows that this quantity 
relaxes to its equilibrium value. 


1 Introduction 

In the last two decades, the problem of dissipation in quantum mechanics, i.e. the consistent 
description of open quantum systems, was investigated by various authors [1, 2, 3, 4, 5]. Because 
dissipative processes imply irreversibility and, therefore, a preferred direction in time, ’f is gen- 
erally thought that quantum dynamical semigroups are the basic tools to introduce dissipation 
in quantum mechanics. In the Markov approximation the most general fornl of the generators 
of such semigroups was given by Lindblad [6]. This formalism has been studied for the case of 
damped harmonic oscillators [7, 8, 9] and applied to various physical phenomena, for instance, the 
damping of collective mddes in deep inelastic collisions in nuclear physics [10] ajnd the interaction 
of a two-level atom with the electromagnetic field [11]. 

In the present work, also dealing with the damping of the harmonic oscillator within the Lind- 
blad theory for open quantum systems, we will explore the physical aspects of the Fokker-Planck 
equation which is the c-number equivalent equation to the master equation for the density opera- 
tor. Generally the master equation gains considerably in clarity if it is represented in terms of the 
Wigner distribution function which satisfies the Fokker-Planck equation. It is worth mentioning 
that these master and Fokker-Planck equations agree in form with the corresponding equations 
formulated in quantum optics [12, 13, 14, 15, 16]. 
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The content of the paper is arranged as follows. In Sec. 2 we review the derivation of the 
master equation of the harmonic oscillator. In Sec. 3 we transform the master equation into the 
Fokker-Planck equation by means of the well-known methods [17, 18, 19]. Then the Fokker-Planck 
equation for the Wigner distribution, subject to either the Gaussian type or the 6-function type 
of initial conditions, is solved by the Wang-Uhlenbeck method. Sec. 4 derives an explicit form of 
the density operator involved in the Lindblad master equation, formulates the entropy using the 
explicit form of the density operator and discusses its temporal behaviour. Finally, concluding 
remarks are given in Sec. 5. 

2 Master equation for the damped harmonic oscillator 

The rigorous formulation for introducing the dissipation into a quantum mechanical system is that 
of quantum dynamical semigroups [2, 3, 6]. According to the axiomatic theory of Lindblad [6], 
the usual von Neumann-Liouville equation ruling the time evolution of closed quantum systems 
is replaced in the case of open systems by the following equation for the density operator p : 

= !(♦.«)• (i) 

Here, denotes the dynamical semigroup describing the irreversible time evolution of the open 
system in the Schrodinger representation and L the infinitesimal generator of the dynamical semi- 
group <J>*. Using the structural theorem of Lindblad [6] which gives the most general form of the 
bounded, completely dissipative Liouville operator L, we obtain the explicit form of the most 
general time-homogeneous quantum mechanical Markovian master equation: 

~~ = L( P (t)) = -~{H,p(t)\ + v+\ + M,p(t)v;)). ( 2 ) 

j 

Here H is the Hamiltonian of the system and the operators Vj and Vj + are bounded operators on 
the Hilbert space of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature are 
of this form after some rearrangement of terms, even for unbounded Liouville opierators. In this 
connection we assume that the general form of the master equation given by (2) is also valid for 
unbounded Liouville operators. 

In this paper we impose a simple condition to the operators H , Vj, Vj + that they are functions 
of the basic observables q and p of the one-dimensional quantum mechanical system (with [<7,p] = 
ih) of such kind that the obtained model is exactly solvable. A precise version for this last 
condition is that linear spaces spanned by first degree (respectively second degree) noncommutative 
polynomials in q and p are invariant to the action of the completely dissipative mapping L. This 
condition implies [7] that Vj are at most first degree polynomials in q and p and H is at most a 
second degree polynomial in q and p. Then the harmonic oscillator Hamiltonian H is chosen of 
the form 

H = H„ + ^(qp + pij), Ho = (3) 


276 


With these choices the Markovian master equation can be written [8]: 


$ = -j[Ho,p] - ^(A +fi)[q,pp + pp] + ^-(A - n)[p,pq + qp] 


dt 


2 ft 


^[ 9 , W.pll - ^Mp, Ip, pH + Ip, pH + Ip, K.pM), 


( 4 ) 


where D pp ,D qq and D pq are the diffusion coefficients and A the friction constant. They satisfy the 
following fundamental constraints [8]: 


A 2 h 2 


i) D pp > 0, ii) D q< i > 0, iii) D pp D qq — D pq > ^ 

In the particular case when the asymptotic state is a Gibbs state 

Pg{o o ) = e -r ^/Tre _E ^, 


these coefficients reduce to 

A + p 


D pp — 


hui A — fi h hu 

hmu> coth D a a = — coth ntrT ^, D pq - U, 


2 2kT' qq 2 mu 

where T is the temperature of the thermal bath 


2kl " 


(5) 


(6) 


(7) 


3 Wigner distribution function 

One useful way to study the consequences of the master equation (4) for the density operator 
of the one-dimensional damped harmonic oscillator is to transform it into more familiar forms, 
such as the equations forj the c-number quasiprobability distributions Glauber P, antinormal 
ordering Q and Wigner W associated with the density operator [20]. In this case the resulting 
differential equations of the Fokker-Planck type for the distribution functions can be solved by 
standard methods [17, 19, 21] employed in quantum optics and observables directly calculated as 
correlations of these distribution functions. 

The Fokker-Planck equation, obtained from the master equation and satisfied by the Wigner 
distribution function W{x\,x 2 ,t) of real variables x\,x 2 corresponding to the operators q,p 


X\ 



x 2 


\/2 h 


mu 


has the form [20]: 


dW 

dt 






t, 7=1,2 


d 2 

dx{dxj 


W, 


where 




1_ / mu)D qq D pq \ 
h V D pq Dpp/ mu ) ' 


( 8 ) 

(9) 


( 10 ) 
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Since the drift coefficients are linear in the variables x x and x 2 and the diffusion coefficients are 
constant with respect to x\ and x 2 , Eq. (9) describes an Ornstein-Uhlenbeck process [22, 23]. 
Following the method developed by Wang and Uhlenbeck [23], we shall solve this Fokker-Planck 
equation, subject to either the wave-packet type or the (5-function type of initial conditions. 

1) When the Fokker-Planck equation is subject to a Gaussian (wave-packet) type of the initial 
condition ( x\q and x 2 o are the initial values of x\ and x 2 at t — 0, respectively) 

W w (x x ,x 2 , 0) = T exp{— 2[(aq - x XG f + (x 2 - x 20 ) 2 ]}, (11) 

the solution is found to be 


W w (x u x 2 , t ) = p=exp{- ^-[<j> w (xi - xj) 2 + i) w {x 2 - x 2 f + Xw{%\ - *i)(x 2 - x 2 )]},(12) 

where 

B w =gig 2 -^gl, g t = g* = ^e 2A< + ^-(e 2At - 1), g 3 = 2[e~ 2Xt + y(l - e“ 2A ')], (13) 

9\ a ■+■ g 2 a — 0«/ — g\ + g 2 ~~ gz, Xw — ^(gi°* 4- g 2 °) ~~ g^i® 4* o*). ( 14) 

We have put a = (g, — iSl)/u>,A = —A — iSl and d x = ( a 2 mu>D qq -f 2 aD pq + D pp / mu) /ft , d 2 = 

( mujD qq + Dp q /to + D pp / rruv)/h and ft 2 = u; 2 — // 2 . The functions x\ and x 2 ^ which are also 

oscillating functions, are given by 

Xj = e _A ‘[x 10 (cosfif + ^ sinOt) + x 20 ^sinfit], (15) 

= e _A< [x 20 (cos Sit — ^ sin Sit ) — -^lo^ sin Sit]. (16) 

2) If the Fokker-Planck equation (9) is subject to the (5-function type of initial condition, the 
Wigner distribution function is given by 

W(x l ,x 2 ,t) = —j—j= exp {--p[M x \ ~ *i) 2 4 tpd(x 2 - x 2 ) 2 + * d (xi - Xi)(x 2 - x 2 )]}, (17) 

where 


B-hh-fl /,=/,• = |(e MI - i), /„=£( 1 (18) 

4>d = fia* 2 + f 2 a 2 - 2/ 3 , ^ = /,+/ 2 - 2 / 3 , Xd = 2[f ia * + f 2 a - f 3 (a + a*)]. (19) 

So, the Wigner functions are 2-dimensional Gaussian distributions with the average values x } and 
x 2 and different widths. 

When time t — > 00 , Xi and x 2 vanish and we obtain the steady state solution: 


W{Z " X2) = 2» N /deL»'(oo) eXpl ~^ i ( ° ll ' ) " l( °° >X ' X ’ 1 - 

The stationary covariance matrix < 7 ^( 00 ) can be determined from the algebraic equation 

A<j w (oo) + <t w (oo)A t = Q W . 


( 20 ) 


( 21 ) 
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4 Entropy and effective temperature 

Entropy is a quantity which may be visualized physically as a measure of the lack of knowledge 
of the system. When we denote by p(t) the density operator in the Schrodinger picture for the 
harmonic oscillator, the entropy S{t) is given by 

S(t) = — &Tr(pln p). (22) 

For calculating the entropy we shall compute straightway the expectation value of the logarithmic 
operator < In p >= Tr(p In p). Accordingly, the problem amounts to derive the explicit form of 
the density operator for the damped harmonic oscillator. 

To get the explicit expression for the density operator, we use the relation p = 27rfiN{VF s (<7,p)}, 
where If , is the Wigner distribution function in the form of standard rule of association and N 
is the normal ordering operator [17, 24] which acting on the function W s (q,p) moves all p to 
the right of the q. By the standard rule of association is meant the correspondence p m q n — ► 
q n p m between functions of two classical variables (q,p) and functions of two quantum mechanical 
canonical operators (q,p). ]The calculation of the density operator is then reduced to a problem of 
transformation of the Wigder distribution function by the N operator, provided tfjat W s is known. 
A special care is necessary for the N operation when the Wigner function is in the exponential form 
of a second order polynomial of q and p. The Wigner distribution function previously obtained 
corresponds however to the form of the Weyl rule of association [25]. The solution (12) of the 
Fokker-Planck equation (9), subject to the wave-packet type of initial condition (11) can be written 
in terms of the coordinate and momentum as: 


W(q,p,t) = — ' l —exp{--^[4>[q- < q >f + < P >f ~ 2 \{q- < q >)(p~ < P >)]}, (23) 

c 2 o 


2tt\/S 


where 


2 h 


< q >= J — — i*| , < p >— v7/mwr 2 , 


(24) 


rmj 


<p = cr p p =< q 2 > - < q > 2 - 


hu 2 ' 1 

rn<jj 


’ V"' W 1 


(25) 


-2 - 2 

ip = < 7 „ =< p l > - < p> = - — 


(26) 


1 /iu> 2 

V = <Tp,(<) = 2 < 9P + P9 > - < « ><P >= gjpX*' i = H>-X 


(27) 


and < A >= Tr(pA) denotes the expectation value of an operator A. The Wigner distribution 
function (23) can be transformed into the form of standard rule of association [26] by 


1 d 2 

W 3 {q,p) = exp(-ik-Q^)W(q,p). 


(28) 
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Upon performing the operation on the right-hand side, we get the Wigner distribution function 
W a , which has the same form as the original W multiplied by h but with x — *h/ 2 in place of x- 
The normal ordering operation of the Wigner function W s in Gaussian form can be carried out by 
applying McCoy theorem [24, 27, 28]. The explicit form of the density operator is the following: 


' = ^ exp[ K -W “ 


cosh *(1 + 


2(£ - 


77) 


x {<t>{q- < q >) 2 + 0(p- < p >) 2 - (x' + *x)[2(?- < q >)(p- < p >) - ih]}]. 


(29) 


where 


£ = <fnl>- x' 2 , x' = X - (30) 

The density operator (29) is in a Gaussian form, as was expected from the initial form of the 
Wigner distribution function. While the density operator is expressed in terms of operators q and 
p, the Wigner distribution is a function of real variables q and p. When time t goes to infinity, the 
density operator approaches to 


p(oo) = 


T 

where a = a pp {oo)cr qq {oo) - cr 2 ? ( oo) and [8]: 

1 


eXp [~2 Ufa ln 2v^- + a ^°°^ 2 ~ + pq)], 


(31) 


<Tpp( ° 0) = 2A(A2+a; 2 -^) ((ma,)V£> ^ + (2A(A "" /i) +u;2 ) Z) pp “ 2rnu 2 (\ - p)D pq ), (32) 


<T, g (°o) 


1 

2(mw) 2 A(A 2 4- u> 2 


-^-((mu;) 2 (2A(A + p) 4- u 2 )D qq 4- u; 2 D pp 4- 2rmo 2 (\ 4- p)D pq ), (33) 
P ) 


<7 P q(oo) 2mA(A 2 4- u> 2 — p 2 )^ (^ “^" ^)( mw ) D qq + (\ p)D pp -\-2m(X p )D pq ). (34) 
In the particular case (7) 

, , h hu hmu> , hu , 

a qq (oo) — 2 mu) 2kT ’ <Tpp ^°°) = 2 co ^ 2kT ’ (7pq ^ <x> ^ = ® (35) 

and the asymptotic state is a Gibbs state (6): 


Pg(oo) = 2 sinh 


Tiuj , 1 1 raw 2 

2TT exp [— +-J " 


«’)]■ 


(36) 
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Because of the presence of the exponential form in the density operator, the construction of the 
logarithmic density is straightforward. In view of the relations (25-27), the expectation value of 
the logarithmic density becomes 


1 h 2 

< In p >= \nh — - ln(6 — — ) 


y/6 2\Jo h 
ti n 2y/a — k 


(37) 


By putting hu = y/S - h/2, we finally get the entropy in a closed form: 

S(t) - k[(v + l)ln(tz + 1) - i/lni/]. (38) 


It is worth noting that the entropy depends only upon the variance of the Wigner distribution. 
When time t — ► oo, the function v goes to s = u)(d%/\ 2 — |di| 2 /(A 2 + fl 2 )) 1/,2 /2ff — 1/2 and the 
entropy relaxes to its equilibrium value 5 ( 00 ) = fc[(s-f 1) ln(s + 1) — s Ins]. It should also be noted 
that the expression (38) has the same form as the entropy of a system of harmonic oscillators 
in thermal equilibrium. In the later case u represents, of course, the average of the number 
operator [29]. While the formal expression (38) for the entropy has a well-known appearance, the 
form of the function v displays clearly a specific feature of the present entropy. We see that the 
time dependence of the entropy is represented by the damping factor exp(— 2\t) and also by the 
oscillating function sm 2 (Qt). The entropy relaxes to its equilibrium value S( 00 ). 


5 Concluding remarks 

Recently we assist to a revival of interest in quantum Brownian motion as a paradigm of quantum 
open systems. There are many motivations. The possibility of preparing systems in macroscopic 
quantum states led to the problems of dissipation in tunneling and of loss of quantum coher- 
ence (decoherence). These problems are intimately related to the issue of quantum- to-classical 
transition. All of them point the necessity of a better understanding of open quantum systems 
and all requires the extension of the model of quantum Brownian motion. The Lindblad the- 
ory provides a selfconsistent treatment of damping as a possible extension of quantum mechanics 
to open systems. In the present paper we have studied the one-dimensional harmonic oscillator 
with dissipation within the framework of this theory. From the master equation of the damped 
quantum oscillator we have derived the corresponding Fokker-Planck equation in the Wigner IV 
representation. The obtained equation describes an Ornstein-Uhlenbeck process. By using the 
Wang-Uhlenbeck method we have solved this equation for the Wigner function, subject to either 
the Gaussian type or the ^-function type of initial conditions and showed that the Wigner func- 
tions are two-dimensional Gaussians with different widths. Then we have obtained the density 
operator. The density operator in a Gaussian form is a function of q, p in addition to several 
time dependent factors. The explicit form of the density operator has been subsequently used to 
calculate the entropy. It relaxes to its equilibrium value. 
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